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Cosmic Microwave Radiation Anisotropics in Brane Worlds 



m 
o 
o 

O 

Q 



> ■ 

oo : 
o . 

m ■ 

o : 

co . 

o ■ 

6 ■ 
a: 



13 



Kazuya Koyama 

Department of Physics, University of Tokyo 7-3-1 Hongo, Bunkyo, Tokyo 113-0033, Japan 

(Dated: February 2, 2008) 

We propose a new formulation to calculate the Cosmic Microwave Background (CMB) spectrum 
in the Randall Sundrum two-branes model based on recent progresses in solving the bulk geometry 
using a low energy approximation. The evolution of the anisotropic stress imprinted on the brane by 
the 5D Weyl tensor is calculated. An impact of the dark radiation perturbation on CMB spectrum 
is investigated in a simple model assuming an initially scale-invariant adiabatic perturbations. The 
dark radiation perturbation induces isocurvature perturbations, but the resultant spectrum can be 
quite different from the prediction of simple mixtures of adiabatic and isocurvature perturbations 
due to Weyl anisotropic stress. 
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1. Introduction 

By suggestions from string theory/M-theory, much at- 
tention is paid to brane world ideas where we are liv- 
ing on a 3-brane in higher-dimensional spacetime. In 
these brane world models, only gravity propagates in 
a higher-dimensional "bulk" spacetime while standard 
model fields are confined to the brane. The simplest re- 
alization of this idea was given by Randall and Sundrum 
lj. The 5D bulk has a negative cosmological constant 
A5 = —6/1 where I is the curvature radius in the bulk 
and the brane has a tension A = 6/k|/ where k| = 8irG5 
and G5 is the 5D Planck constant. 

In this model, predictions of 4D general relativity are 
modified due to the influence of the gravitational fields 
in the bulk. These modifications need to be consis- 
tent with cosmological observations. Among them, the 
observations on Cosmic Microwave Background (CMB) 
anisotropics are now dramatically improving. The 
WMAP experiment has alreadyprovided precise mea- 
surements of CMB anisotropies |2|. It is then necessary 
to calculate the CMB spectrum in brane worlds in order 
to check the consistency of the model. 

These works were initiated soon after the paper of Ran- 
dall Sundrum |3j. Although a large number of papers 
have tried to give a prediction of CMB anisotropy, there 
is still no quantitative estimation. This is due to the 
difficulties in solving full 5D perturbations Q. In view 
of observational improvements, it is eagerly desired to 
develop a formulation which enables us to provide quan- 
titative predictions. In this Letter, we propose such a 
formulation based on recent progresses in solving bulk 
geometry using a low energy approximation 

We consider a simple Randall-Sundrum type two- 
branes model where a low energy approximation is ap- 
plicable. The stabilization mechanism is not introduced 
and the physical brane is assumed to be the positive ten- 
sion brane. This simple setup enables us to understand 
clearly how the bulk gravitational fields affect the evolu- 
tion of perturbations on the brane. 



2. View from the brane 

The effective 4D Einstein equation on the brane is given 
by@: 
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where k\ = 8irG4 = and is the quadratic 

function of the energy momentum tensor, which can 
be neglected at low energies Tjf/\ <C 1. The 
energy-momentum tensor satisfies the conservation's law 
V^T^jy = where V M is the covariant derivative on the 
brane. E^ is the electric part of the 5D Weyl tensor and 



carries the information in the bulk, 
identities and the symmetry of Weyl tensor, 



satisfy 



= 0, 



From the 4D Bianci 
£,,„ should 
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on the brane for T M „/ 'A <C 1. An important point is that 
£ M „ cannot be determined only from equations on the 
brane Eq. ■ One needs 5D equations for £^ v in the bulk 
and this is the source of great complexity of the problem. 
To observe this fact, it is convenient to parameterize E^ v 
as an effective energy-momentum tensor; 
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where Y(k,x) oc e lkx is the normalized scalar harmonics 
and the vector Yi and traceless tensor Yja are constructed 



from Y as K- 
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<-lj — n T SijY/3. 

In the background universe, Eq.<0 gives p £ +4Hp £ = 
and P £ — p £ /3, where H = a/a and a is the scale factor. 
The solution for Weyl energy density is then given by 
p £ = Ca~ 4 where C is the integration constant. Thus 
the contribution from E^ u appears as dark radiation. It 
has been shown that C is related to the mass of the AdS- 
Schwarzshild BH in the bulk. In the following, we will 
assume p £ — in the background spacetime. 
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Let us consider the perturbations around the above 
background spacetime. The linear scalar metric is taken 
as 

ds 2 = -(1 + 2^{t)Y)dt 2 + a{tf{\ + 2<ff{t)Y)S lj dx l dx j . 



Eq.© give 
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-a- x W £ , SP £ = -8pe, 
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for p/X <C 1. At large scales ka~ 1 /H — > 0, we can neglect 
the terms proportional to k. The solution for Sp £ is given 
by Sp £ = 8Ca~ A where 5C is the integration constant. 
The amplitude SC is related to the perturbatively small 
AdS-Schwarzshild mass in the bulk. A problem is that it 
is impossible to determine anisotropic stress Stt £ because 
it is dropped from Eq.Q for kaT 1 jB — ► 0. In |lCj. it 
was clearly shown that this uncertainty of Stt £ prevents 
us from predicting CMB anisotropy. The temperature 
anisotropy caused by (direct) Sachs- Wolfe effect is given 
by the formula (lo| 
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where Oo is the temperature anisotropy of radiation and 
we used Eq.QJ in the second equality. The evolution of 
C = $ + 5p/3(p + P) is determined only by the conser- 
vation of the matter energy-momentum tensor. We will 
assume the matter perturbation is adiabatic SP = c 2 Sp 
where Sp is the density perturbation, SP is the pressure 
perturbation and c s is the sound velocity. The continuity 
equation for matter perturbation implies £ = = const. 
In addition to £, the solutions for metric perturbations 
$ and 4" are needed. From the effective 4D equation Q), 
we can get the equations to determine <I> and iff as 



Got = + * = -n\k- 2 a 2 ^ £ , (5) 



where £ tot is the curvature perturbation on hypersurface 
of uniform total energy density; 



Ctot — c 
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Note that the solution for C, to t was obtained only by 
the law of conservation. However, the latter equation 
in Eq.|0 contains undetermined 5ir £ . Indeed, Eq.|@} is 
written at the decoupling time as 

AT I 2 ( P r\, r 

- nlk~ 2 a 2 5TT £ + 2 K 2 4 k~ 2 a- 5/2 J a 7/2 STT £ da,(7) 



where <5C* = Sp £ / p r and p r is the radiation energy den- 
sity. The second term comes from the entropy perturba- 
tion contribution due to Weyl energy density. Unless the 
behavior of 5tt £ is known, we cannot say anything about 
the effect of this dark radiation. It is possible to assume 
some ansatz for 5tt £ , say Stt £ oc Sp £ However, be- 

cause 5tt £ is induced by 5D gravitational perturbations, 
it is essential to solve the evolution equations for £ M „ in 
the bulk. 

3. Solution for 8^ v in the bulk 

The evolution equations for obtained from 5D 

Bianchi identity are consisted of discouragingly compli- 
cated partial differential equations. Fortunately, the cur- 
vature radius on the brane defined by L 2 ~ 1/V„ ~ 
1 / n\Tff at the decoupling time is significantly longer than 
the curvature radius in the bulk I because the experi- 
ments on Newton's force law already impose the con- 
straint on I as / < 1 mm. Thus the system has a natural 
small parameter e = (l/L) 2 ~ T£/\ and it is possible 
to solve the equations by the perturbation in terms of 
e. If we consider a second brane with negative tension 
A c = — 6/n 2 l at the physical distance dg(x)l from our 
brane, £^ v can be determined as 



£,? = 
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where T£ v is the energy-momentum tensor on the second 
brane. Because the bulk spacetime shrinks exponentially 
due to the negative cosmological constant in the bulk, 
the curvature radius on the second brane L c is shorter 
for larger do as L c — Le~ d<) . Thus in order to ensure 
that the low energy approximation can be applied on the 
second brane (l/L c ) 2 < 1, the radion do should satisfy 
e- 2do > {l/L) 2 < 1. 

The evolution of do can be calculated from the traceless 
condition ££ = in Eq.0; 



V 2 d -(Vd ) 2 = 5(T + e- 2d «T c 
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(9) 



From Eqs.(|SJ and the behavior of £(} is completely de- 
termined. Substituting the solution for into Eq. 0J, 
the effective theory on the brane becomes quasi-scalar- 
tensor theory [l^ . 

4- View from the bulk 

The solution for (Eq.JSJ)) should be consistent with 
Eq.©. In the background spacetime, p £ is written in 
terms of do and energy densities on both branes; 
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The evolution equation for do can be integrated once with 
the help of the energy-momentum conservation on each 
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brane. We get p £ — Ca~ A where C is the integration con- 
stant. In this case, the integration constant is interpreted 
as the initial condition for the radion do(t). Now we turn 
to the perturbations. We denote the perturbation of the 
radion as do(t, x) = do(t)+N(t)Y(k, x). Sp £ can be eval- 
uated in the same way as the background spacetime. We 
get 5pe = 5Ca~ 4 at large scales where SC is the integra- 
tion constant associated with N. An advantage of our 
approach is that an equation for Weyl anisotropic stress 
Sir £ can be derived as 



gives the evolution equation for E; 
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Hence the behavior of Stt £ is determined by the radion 
perturbation TV. Because the radion perturbation N is 
coupled to the metric perturbations $ and \P, we should 
solve them at the same time. The behavior of the metric 
perturbations is determined by imposing the equation of 
state on the matter perturbations: SP — c 2 5p and 5P C = 
c 2 c Sp c . We then have three equations, SP = c 2 Sp, SP C = 
c 2 c Sp c and $ + \P = — K 2 k~ 2 a 2 STr £ for three unknown 
functions N, $ and $. It is convenient to introduce a 
new set of variables tp, ip c and E to solve these equations 
@. We take 



The function E is identified with the bulk anisotropic 
perturbation We have a closed set of equations for 
ip, ip c and E. Therefore, the solution for Sirs can be ob- 
tained from Eas. flU)! and JTJJ. 

5. CMB anisotropy in a simple model 
Let us consider the simplest case where do = d* = const. 
and d* is sufficiently large e~ 2rf * <C 1. It can be real- 
ized by taking p c = —pe 2d * and w = w c . This choice 
is consistent with p £ = 0. First, let us consider large 
scale perturbations. Assuming the scale factor is given 
by a oc t 2 / 3 ( 1 +"') ( w =const.), the solution for ip is ob- 
tained as 
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E, N — <p c - ip ~ k~ 2 a 2 d Q E, 
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The equation SP = c 2 Sp gives the equation for ip; 



p+{2 + 3c 2 s )Hp - (3H 2 + 2H + 3c 2 H 2 )ip 



for w ^ —1/3, 1/9, where £ c * is the curvature perturba- 
tion on hypersurface of uniform matter energy density on 
the second brane. Here we neglected the homogeneous 
solutions. Note that w = 1/9 is not singular and one can 
find a solution for ip as — {l/2)5C*{p/ p r ) lna. One point 
is that ip c that depends on SC* satisfies <p c = e 2d *ip, then 
does not contribute to E. 



2c SC 

4 solution for metric perturbations 



ip- 



E 



(12) 



Now we can obtain the 
We find the parts of $ 
and \& corresponding to each term in the solution for ip: 



At large scales ka~ 1 /H — > 0, the equation for ip is de- 
coupled from E and it can be integrated once to give the 
first order differential equation for ip; 
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Here, and are the same as conventional 4D solu- 
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where C* is the integration constant. We should note 
that the left-hand side of Eq. l|13|) is nothing but the 
solution for ( tot (see Eqs.ljHJ) and CU)- Then the be- 
havior of ip determines the evolution of the curvature 
perturbation which is independent of do- The evolu- 
tion equation for p c is obtained from SP C — c 2 c Sp c 
and it is given by replacing H,a,c 2 ,d/dt, and Sp £ to 
H c = e d "(H - d ),a c = ae~ d ° , c 2 sc , e d °d/dt and e Ad °5p £ , 
respectively, in Ea. l|12(l . The function ip and ip c describe 
the displacement of our brane and the second brane re- 
spectively, and their relative difference causes the radion 



perturbation N. The equation $ + 



-K%k 2 a Sirs 



tions except for the additional terms which depend on 
the radion d* and "shadow matter" [||. These addi- 
tional terms can be neglected for e~ 2d * <C 1. Then we 
find the solution for Weyl anisotropic stress as 

nlk- 2 a 2 8TT £ = 2tp £ . (16) 

Large scale CMB anisotropy can be written as 

Dark radiation does not affect AT/T. This is quite a 
non-trivial result. As mentioned before, the Weyl energy 
density perturbations induce the entropy perturbations. 
The effect of Weyl anisotropic stress Sir £ exactly cancels 
this entropy perturbation in Eq.Q. 
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Let us investigate whether this cancellation holds for 
small scale perturbations. Under a tight coupling approx- 
imation of baryon-photon fluid, the evolution equation 
for radiation temperature anisotropy Qq becomes |l3f 



HQr 



R 



1 + R 



Q + k 2 a~ 2 c 2 s Qo = F, 



R 



F = -HQ- - y«- 2 f, (18) 

where R = 3pb/4:p r and pb is the baryon energy den- 
sity. The solutions for metric perturbations are given by 
the upper equation in Eq. 1)1 5JI though ips is the solutions 
for Ea. (|12fl with E = and the equation for Sps is ob- 
tained by Eqs.© and l|l(j|) . The source term F can be 
decomposed as F = F^ + Fg according to the solutions 
for metric perturbations (Eq. i|15|) ). Observed tempera- 
ture anisotropy (Eq.lgJ) is given by AT/T = 6 + * = 
©o — f£ + Thus if Q 0e = ©o — <Pe behaves in the 
same way as conventional 4D theory, the temperature 
anisotropy also behaves in the same way. From Ea. (ll8ll . 
we get the equation for <3q £ : 



e 0e 



HQ 



R 



0e ' 



1 + R 



Q ne + k 2 a 2 c 2 Q 0e 



1 



(19) 



Hence at large scales k — > or in the radiation domi- 
nated universe c 2 — 1/3, the temperature anisotropy is 
exactly the same as the conventional 4D theory. But as 
the Universe becomes matter dominated c 2 ^ 1/3 and 
perturbations enter the horizon, the behavior of temper- 
ature anisotropy is modified. Because the amplitude of 
(pe decreases with time under the horizon, the effect is 
largest for the first acoustic oscillation. Fig.l shows the 
resultant CMB anisotropy for given cosmological parame- 
ters. One can see that the dark radiation does not modify 
the CMB anisotropy at I ~ 2 and at large I but it mod- 
ifies a first acoustic oscillation. The resultant spectrum 
is quite different from the prediction of simple mixtures 
of adiabatic and isocurvature perturbations due to Weyl 
anisotropic stress 5tt£. 

6. Discussions 

In this letter, we developed a formulation to calculate 
CMB anisotropy in two branes model at low energies with 
an appropriately small distance between two branes. In 
realistic models, we should introduce stabilization mech- 
anisms of the radion. Our formulation can be extended 
to include stabilization mechanisms. We considered only 
the classical theory of the perturbations. Thus the ini- 
tial spectrum remains to be determined, and it can be 
modified by the effects from the bulk in the brane world. 
These issues need further investigations. 



[1] L. Randall and R. Sundrum, Phys. Rev. Lett. 83, 3370 
(1999); ibid, 83, 4690 (1999). 



[2] 
[3] 



[1] 



[5] 

[6] 
[7] 



[9] 
[10] 

[11] 
[12] 



[13] 
[14] 



C.L. Bennett et. al, astro-ph/0302207 

S. Mukohyama, Phys. Rev. D 62, 084015 (2000); H. Ko- 

dama, A. Ishibashi and O. Seto, Phys. Rev. D 62, 064022 

(2000) ; R. Maartens, Phys. Rev. D 62, 084023 (2000) 

; D. Langlois, Phys. Rev. D 62, 126012 (2000) ; C. van 

de Bruck, M. Dorca, R. H. Brandenberger and A. Lukas, 

Phys. Rev. D 62, 123515 (2000) ; K. Koyama and J. 

Soda, Phys. Rev. D 62, 123502 (2000) . 

see, for example, N. Deruelle, Proceedings of the Porto 

JENAM conference (gr-qc/0301035 ); R. Maartens, Prog. 

Theor. Phys. Suppl. 148, 213 (2003). 

C. Csaki, M. Graesser, L. Randall and J Terning, Phys. 

Rev. D62, 045015 (2000); 

K. Koyama, Phys. Rev. D 66 084003 (2002). 

T. Wiseman, Class. Quantum Grav. 19, 3083 (2002); S. 

Kanno and J. Soda, Phys. Rev. D66,083506 (2002). 

T. Shiromizu and K. Koyama, Phys. Rev. D67, 084022 

(2003). 

T. Shiromizu, K.I. Maeda and M. Sasaki, Phys. Rev. 
D62, 024012 (2000); 

D Langlois, R Maartens, M Sasaki, D Wands, Phys. Rev. 
D63, 084009 (2001). 

J. D. Barrow and R. Maartens, Phys. Lett. B532, 153 
(2002). 

The CMB anisotopies in scalar-tensor theory have been 
investigated intensively, see for example, R. Nagata, 
T. Chiba and N. Sugiyama, Phys. Rev. D66, 103510, 
(2002); X. Chen and M. Kamionkowski, Phys. Rev. 
D60, 104036 (1999). However, the initial conditions for 
scalar field perturbations in these literatures correspond 
to <5C* = in the brane world context. Thus we cannot 
apply their results for the investigation of the effect of 
dark radiation perturbation. 

W. Hu and N. Sugiyama, Astrophys. J. 444, 489 (1995). 
M. Doran, astro-ph/0302138 U. Seljak and M. Zaldar- 
riaga, Astrophys. J. 459, 437 (1996). 





A 


Sp E /pr =5C- 


m 


5C.= -2 Cf 

5C-= -£. 

8C.= 

5C.= £. 

8(7.= 2 Cf 

j WMAP 


a\ 

iijij \ 

Jj w 















FIG. 1: CMB angular power spectrum for various SC. £* 
is appropriately normalized. We take Qq = 0.26, S1a = 
0.70, Q b = 0.04, h = 0.72 and n = 1. The observational data 
was taken from WMAP Q and the spectrum was calculated 
by modifying CMBEASY which is based on CMBFAST Q. 



